ARID 


doi; 10. 13866/;. azr. 2019.02.01 


X ive f'ERHT 


ChinaXi HHT 


旱 区 研究 
ZONE RESEARCH 


1971—2014 年 青藏 高 原 参 考 蒸 散 变化 及 其 归 因 
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摘 X: 人 研究 参考 蒸 散 时 空 变化 格局 并 辨识 其 驱动 因子 ,是 认识 区 域 水 文 过 程 及 其 对 气候 变化 响应 的 重要 途径 。 
基于 修正 的 FAO 56 Penman-Monteith 公式 和 青藏 高 原 75 个 台 站 逐日 气象 观测 资料 ,分 析 了 1971 一 2014 年 高 原 参 
考 蒸 散 变 化 的 转折 特征 ,并 探讨 转折 前 后 的 年 际 与 季节 变化 趋势 及 其 主导 因素 。 结 果 表 明 :1971 一 1996 年 青藏 高 
原 参 考 蒸 散 呈 急剧 下 降 态 势 [ - 27. 07 mm + (10a) ~), fii 1997—2014 年 高 原 参 考 敬 散 增加 趋势 显著 [40. 16 
mm * (10a) ^) , 尤 以 33"N 以 南 区 域 最 为 突出 。 这 与 影响 参考 燕 散 变化 的 气候 因子 变化 趋势 的 年 际 转折 密切 相 
关 。 其 中 ,风速 变 小 是 1971—1996 年 高 原 年 参考 蒸 散 减 少 的 主要 因素 ,特别 是 在 高 原 北 部 ;相对 湿度 降低 则 极 大 
地 促进 了 1997—2014 年 高 原 主体 ( 除 高 原 北 缘 外 ) 参 考 蒸 散 的 显著 增加 。 此 外 ,从 季节 上 看 ,春秋 、 冬 季 参 考 燕 散 
变化 的 最 大 贡献 因子 由 之 前 的 风速 减 小 转变 为 1997—2014 年 的 相对 湿度 下 降 ;影响 高 原 夏 季 参 考 燕 散 的 主导 因 
子 是 1971 一 1996 年 相对 湿度 的 增加 ,之 后 则 转变 为 1997 一 2014 年 日 照 时 数 的 增加 。 
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气候 持续 变 暖 导致 全 球 水 循环 加 剧 "] 。 作 为 
水 循环 过 程 的 重要 环节 之 一 , 蒸 散 直接 决定 了 地 表 
和 大 气 间 的 水 汽 传输 。 人 参考 蒸 散 表 示 在 一 定 气 象 条 
件 下 水 分 供应 不 受 限 制 时 某 一 固定 下 热 面 可 能 达到 
的 最 大 莹 散 量 ,这 为 评价 和 比较 不 同 气 候 条 件 下 的 
蒸 散 差 异 提 供 了 一 个 标准 。 参 考 蒸 散 既 是 陆 面 水 分 
平衡 和 地 表 能 量 平衡 的 主要 组 成 部 分 ,又 是 估算 实 
RZE .评价 区 域 干 湿 状况 和 制定 水 资源 管理 计划 
的 关键 因子 2] 。 因 而 ,辨识 全 球 变 暧 背景 下 影响 参 
考 燕 散 变 化 的 的 主导 气候 因子 ,成 为 区 域 水 文 过 程 
研究 的 重要 内 容 , 从 而 有 效 指导 区 域 水 土 配置 。 

青藏 高 原平 均 海拔 超过 4 000 m, 总 面积 约 
2.5.x 10° kn ,是 我 国 重要 的 生态 安全 屏障 。 
其 感 热 加 热 和 潜 热 释放 能 加 剧 季 风 环 流 ,进而 影响 
区 域 气候 系统 5] 。 在 青藏 高 原 气温 加 速 升 高 ““”、 
风速 减 小 5 太阳 辐射 下 降 ” 呈 等 气候 变化 背景 
下 ,高 原 参 考 蒸 散 变化 及 其 归 因 研究 越发 受到 关注 。 
参考 燕 散 变化 归 因 的 研究 方法 包括 相关 分 析 、 去 趋 
3A .导数 敏感 性 系数 等 。 其 中 ,相关 分 析 计 算 燕 发 
和 气候 因子 的 相关 系数 ,或 者 求 取 二 者 的 回归 方 
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FEO” 。 相 关系 数 越 高 ,表明 气候 因子 对 蒸 散 变化 影 
响 越 大 。 去 趋势 方法 通过 评估 气候 因子 的 变化 趋势 
来 甄选 主导 因子 ,去 趋势 后 导致 蒸发 出 现 最 大 差 值 
的 因子 被 视 为 影响 蒸 散 的 主要 原因 与 -  。 基 于 导 
数 方法 , 求 得 茸 散 各 气候 因子 的 导数 ,定量 化 评价 它 
们 对 莹 散 变化 的 贡献 大 小 -9 。 敏 感性 分 析 研 究 
气候 因子 变化 引起 的 蒸 散 模 型 输出 响应 ,通过 假定 
其 他 因子 不 变 , 分 析 单 个 因子 的 变化 对 模型 的 作用 ， 
进而 确定 敏感 性 系数 '""。 目 前 , 归 因 方法 已 被 广泛 
用 于 分 析 青 藏 高 原 潜 在 蒸 散 的 影响 因子 研究 。 
Chen 等 5 研究 认为 ,1961 一 2000 年 促使 高 原 潜 在 
茸 散 减少 的 主要 因子 是 风速 和 相对 湿度 ; Zhang 
457 4 iH ,1971—2004 年 风速 .日照 时 数 减 小 分 别 
是 青藏 高 原 北部 、 东 南部 参考 蔚 散 显著 下 降 的 主导 
因子 ; 尹 云 鹤 等 5 认为 ,风速 减 小 是 影响 青藏 高 原 
年 潜在 蒸 散 降低 的 主要 原因 ;Xie 和 Zhu 指出, 风 
速 和 实际 水 汽 压 是 该 地 区 参考 蒸 散 下 降 的 主要 原 
因 ; 黄 会 平等 通过 主 成 分 分 析 发 现 , 包 括 气温 和 
水 汽 压 在 内 的 热力 学 因子 影响 参考 蒸 散 ,在 青藏 高 
原 地 区 呈 环 形 向 中 心 递 增 ;Zhang 等 3! 则 明晰 了 青 
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藏 高 原 农 业 区 参考 蒸 散 的 敏感 因子 为 实际 水 汽 压 。 
值得 注意 的 是 ,数值 天 气 模式 不 同 于 之 前 的 气象 站 
点 资料 , 其 数据 开始 应 用 于 高 原 参 考 蒸 散 研 
和 

lA ERRE 多 利用 FAO 推荐 的 Penman- 
Monteith 公式 分 析 参 考 蒸 散 时 空 分 布 特征 ,同时 采 
用 偏 相 关 分 析 或 敏感 性 系数 探讨 不 同 气象 要 素 对 参 
考 蒸 散 变化 的 影响 。 但 上 述 研 究 多 选取 整个 研究 时 
段 进 行 趋势 分 析 ,考虑 到 高 原 气 候 变 化 的 突变 特征 ， 
有 必要 检测 研究 时 期 内 的 参考 蒸 散 突变 ,进而 分 析 
突变 前 后 不 同 阶 段 参考 蒸 散 变 化 的 趋势 差异 及 其 原 
因 。 本 文 基于 1971—2014 年 青藏 高 原 75 个 气象 站 
点 的 日 值 资料 ,系统 分 析 过 去 44 a, 特 别 是 突变 前 后 
阶段 的 参考 蒸 散 时 空 变化 特征 ,重点 探讨 参考 蒸 散 变 
化 的 影响 因子 及 其 区 域 分 异 规律 ,为 深入 揭示 全 球 变 
暧 背 景 下 青藏 高 原 区 域 水 文 过 程 提供 科学 依据 。 


1 数据 与 方法 


1.1 气象 数据 来 源 

气象 数据 来 源 于 中 国 气象 数据 网 (http://data. 
cma. cn/ ) ,包括 1971—2014 年 青藏 高 原 及 其 周边 
地 区 75 个 气象 站 点 的 日 平均 气温 、 最 高 气温 、 最 低 
气温 .平均 相对 湿度 、 降 水 量 . 平均 风 速 、 气压 以 及 
16 个 辐射 站 点 的 数据 (图 1) 。 其 中 ,75 个 气象 站 点 
分 布 在 西藏 (22 个 ) 青海 省 (23 个 ) .新 疆 (5 个 )、 
甘肃 省 (8 个 ) 四川 省 (14 个 ) 云南 省 (3 个 ) E 


80°E 


+ 辐射 站 点 
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述 站 点 平均 海拔 高 度 为 3 148.2 m, 海 拔 最 高 .最 低 
的 站 点 分 别 为 安 多 站 (4 800. 0 m) . 若 芜 站 (887.7 
m) 。 考 虑 到 1970s 前 的 站 点 观测 数据 完整 性 较 差 ， 
本 文选 取 1971 一 2014 年 的 资料 。 
1.2 研究 方法 
1.2.1 FAO Penman-Monteith 方法 KWR 
散 模 型 , 使 用 了 1998 年 FAO 提出 的 Penman- 
Monteith 最 新 修正 公式 , 即 FAO Penman-Monteith 模 
型 中 。 计 算 公式 如 下 : 

900 


0. 408A(R, - C) +Y p.273. 一 es) 


More Asyl +0. 34u;) d 


式 中 :RR, 为 净 辐 射 (MJ © m””. d”) ;6 为 土壤 热 通 
量 (MJ +m? + d7)3;7 为 日 平均 气温 (% ) iu, X2 
m 高 度 处 的 风速 Cm * s 7) ;e, 为 他 和 水 汽 压 (kPa ) ; 
e, 为 实际 水 汽 压 (kPa) ;A 为 饱和 水 汽 压 对 气温 了 的 
斜率 (kPa . Co") ;y 为 湿度 计 常 数 (kPa . To") 

(1) 太阳 辐射 校正 

太阳 辐射 通常 采用 Angstrom 公式 进行 校正 ,其 
公式 如 下 中 


R, = (a, +b, 二)R (2) 


式 中 :及 为 太阳 辐射 (MJ m? + d7');R, HRA 
辐射 (MJ . m^ *d^);a, 为 阴 天 天 气 (n 20) KA 
顶 辐 射 比值 ,a, +b, 为 晴天 天 气 大 气 顶 辐射 比值 ;n 


为 每 日 日 照 时 数 ;N 为 最 大 可 能 日 照 时 数 ;六 为 相对 


90*E 100°E 


A1 青藏 高 原 气 象 站 点 分 布 示 意图 


Fig.1 Distribution of meteorological stations over the Tibetan Plateau 
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汪 步 惟 等 :1971 一 2014 4F 6 HR ES 5. 75 RB hI ES] 


日 照 时 数 。 

本 文 利 用 青藏 高 原 地 区 已 有 16 个 辐射 站 点 观 
测 资料 进行 回归 分 析 , 来 确定 系数 a, 和 565., 并 通过 
ArcGIS 9.3 克 里 金 (Kriging) 方 法 对 整个 青藏 高 原 地 
区 进行 系数 插值 ,最 后 提取 75 个 气象 观测 站 点 的 a, 
Fil b, 系数 。 

(2) 净 长 波 辐射 

ADMIT Astin FO : 


4 4 
T max ,上 + T min,k 


5 (0.34 - Je, ) 
(党 一-0.35]] (3) 


式 中 :co 为 斯 提 芬 - 玻 尔 效 曼 常数 (4. 903 x 107 
MJ m”? +d! + K); 7 是 一 天 中 最 大 绝对 温 
JE CK) T, 是 一 天 中 最 小 绝对 温度 (K) ;e, 为 实际 


Ry = | 


ES E 
水 汽 压 (kPa) ;R, 为 大 气 硕 太阳 辐射 ,是 相对 短 
波 辐 射 率 。 
(3) 土壤 热 通 量 
土壤 热 通 量 利用 月 平均 气温 的 经 验方 法 进行 估 
算 HARUT” : 
G nih = 0. 07 ( T ost tl 一 T nih -1 ) (4 ) 
或 者 , Pecos +1 当 未 知 时 , 
G monih,i z 0. 14( T ds = Tatai ) (5 ) 


式 中 :Go 为 第 i 月 的 土壤 热 通 量 (MJ m - 
d^ s T uis osa s asas IRA d - 1 ii + 1 
月 的 平均 气温 。 

(4) 净 辐 射 通 量 

净 辐 射 通 量 的 计算 公式 如 下 中 : 

R,=R,, +R, +G (6) 

SOE LR, AG LG (MJ - m d RUN 
净 长 波 辐 射 通 量 (MJ m” + d7') ;6 为 土壤 热 通 量 
(MJ md)。 
1.2.2 统计 方法 

(1) 趋势 分 析 


ias; det E55) ( E25) 
b= er. (7) 
( xg POSSE 


(2) 滑动 1 检验 

本 文 利用 滑动 上 检验 检测 1971—2014 年 高 原 
BARE .气温 .日照 时 数 、 风速 和 相对 湿度 突变 及 
其 发 生 的 时 间 。 该 检验 是 通过 考察 2 组 样本 平均 值 
的 差 值 是 否 显 著 来 检验 突变 。 其 基本 思想 是 把 
一 气候 序列 中 两 段子 序列 均值 有 无 显著 差异 视 为 来 
自 两 个 总 体 均值 有 无 显著 差异 的 问题 来 检验 。 如 两 
段子 序列 均值 差异 超过 一 定 的 显著 性 水 平 ,可 以 认 
为 均值 发 生 质 变 ,有 突变 发 生 。 对 具有 2 个 样本 量 
的 时 间 序 列 *, 人 为 设置 某 一 时 刻 为 基准 点 ,基准 点 
前 后 两 段子 序列 x 和 x, 的 样本 分 别 为 n, 和 nn ,两 
段子 序列 平均 值 分 别 为 和 名, 方差 分 别 为 51 和 
5;。 定 义 统计 量 为 : 


NL E M (8) 


€ 
. 


n, Ny 


其 中 ,序列 x 的 标准 差 ; 为 : 
nis? *njs? 
in n, tn, 一 2 e 
2 结果 与 分 析 


2.1 青藏 高 原 参 考 蒸 散 变化 趋势 
2.1.1 参考 蒸 散 空间 平均 值 变化 趋势 转变 
1971—2014 年 青藏 高 原 参 考 蒸 散 多 年 平均 值 为 
917. 15 mm , 波动 较 强 ( 图 2a) ,最 大 值 (1 408. 29 
mm)、 最 小 值 (672. 61 mm) 分 别 出 现 在 若 芜 、 错 那 
站 。 

高 原 参 考 蒸 散 累 积 距 平 在 1997 年 前 后 由 正 值 
转 为 负 值 (图 2b) ,结合 高 原 气候 多 因子 突变 分 析 结 
果 ,高原 参 考 蒸 散 变 化 趋势 在 1997 年 前 后 发 生 转 
变 。 据 此 ,本 研究 对 比分 析 了 趋势 突变 前 后 2 个 时 
段 (1971 一 1996 年 和 1997—2014 4E) 高原 参 考 蒸 散 
变化 及 其 驱动 因子 的 区 域 差 异 。1971 一 1996 年 ,高 


为 分 析 1971—2014 年 高 原 各 站 点 参考 蔡 散 \ 气 
温 \ 日 照 时 数 、 风 速 和 相对 湿度 的 线性 变化 趋势 ,本 
文采 用 线性 变化 率 “-” 对 观测 数据 x; 及 相应 的 时 
I t, 进行 处 理 。 线 性 变化 率 " 的 最 小 二 乘法 估算 公 
式 为 : 


原 参 考 蒸 散 减少 趋势 显著 [ -27.07 mm : (10a) ^); 
1997—2014 年 ,高 原 参 考 蒸 散 呈现 显著 的 增加 趋势 
(40.16 mm * (10a) ~!) (I 2a) , 
2.1.2 各 站 点 参考 蒸 散 年 变化 趋势 ”1971 一 2014 
青藏 高 原 各 站 点 参考 蒸 散 年 变化 趋势 空间 分 布 不 
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1971 一 2014 年 青藏 高 原 参考 长 散 趋势 分 析 (a) 和 累积 距 平 (b) 


Fig.2 Trend (a) and cumulative anomaly (b) of annual reference evapotranspiration over the Tibetan Plateau during 


the period of 1971 - 2014 
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(c) 1997—2014 年 


到 3 ”不同 时段 青藏 高 原 各 站 点 参考 蒸 散 线性 变化 趋势 


Fig.3 Linear trends of reference evapotranspiration 
at the meteorological stations over the Tibetan 


Plateau during different periods 


均匀 (图 3a) 75 个 站 点 中 ,年 参考 蒸 散 显著 (P < 
0. 05) 减少 和 增加 的 气象 站 点 分 别 为 15 16 个 ,其 线 
性 变化 率 主 要 介 于 - 49.99 ~ - 10. 00 mm - 
(10a) -10.01 ~50.00 mm * (10a) -之 间 ; 年 参考 
蒸 散 显著 增加 的 站 点 多 集中 分 布 在 巴 颜 喀 拉 山 和 阿 
尼 玛 卿 山 周 边 地 区 。 

1971—1996 年 高 原 参 考 蒸 散 整 体 呈 减少 趋势 ， 
有 31 个 站 点 达到 显著 水 平 (P <0. 05 ) ;其 中 ,线性 
变化 率 < - 100.00 mm . (10a) ^, - 99. 99 ~ 
-50.00 mm + (10a) ', -49.99 ~ -10.00 mm . 
(10a) 的 站 点 个 数 比例 为 4:9:18( 图 3b)。 

相 比 之 下 ,1997 一 2014 年 整个 高 原 参 考 蒸 散 增 
加 趋势 更 为 一 致 :显著 增加 的 站 点 多 达 46 个 , 且 其 
线性 变化 率 普遍 超过 50.01 mm + (10a) ^ ,这 在 青 
藏 高 原 约 33°N 以 南 ( 唐 古 拉 山 以 南 区 域 ) 尤 为 明显 
(图 3c)。 具 体 来 讲 ,显著 增加 的 站 点 中 ,其 线性 变 
化 率 > 100.01 mm * (10a) ^ 50.01 ~100.00 mm - 
(10a) ^! ,10. 01 ~50.00 mm -* (10a) 一 的 站 点 数 比 
例 为 3:24:19。 
2.1.3 各 站 点 参考 蒸 散 季节 变化 趋势 ”突变 前 后 ， 
青藏 高 原 各 季 参 考 蒸 散 变化 趋势 的 区 域 差 异 均 比 较 
显著 (图 4) 。1971 一 1996 年 春 . 夏 、 秋 ,冬季 高 原 75 
个 站 点 参考 蒸 散 的 线性 变化 率 平均 值 为 
-9.18 mm * (1043) !, - 9. 56 mm * (10a)7', 
-4.49 mm * (103) !, -2. 60 mm * (103) ', H 
中 , 春 、 夏 季 参 考 蒸 散 显著 减 小 的 区 域 较为 相似 , E 
要 分 布 在 昆仑 山 一 阿尔 金山 一 祁连山 沿线 和 横断 山 
脉 地 区 ,下 降幅 度 多 在 - 49.99 ~ - 10. 00 mm - 
(10a) -之 间 ( 图 4a.4c) KE ,冬季 参考 蒸 散 整体 
呈 减 少 趋势 ,显著 减少 站 点 个 数 分 别 为 20 .24 ,冬季 
减少 幅度 相对 秋季 较 弱 (图 4e、4g)。 

1997 一 2014 年 春 、 夏 秋 、 冬 季 高 原 75 个 站 点 
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1971 一 1996 年 


25°N 30°N 35N 40N 


80° E 90° E 100° E 


25°N 30°N 35°N 40°N 


80° E 90° E 100° E 


25°N 30°N 35°N 40N 


1997—2014 年 


线性 变化 率 / @ <-100.00 @ -49.99--10.00 90.01-10.00  @50.01~100.00 @ 不 显著 


[mm:.(10a)-] @ -99.99--50.00 O -9.99~0.00 


@ 10.01~50.00  @>100.01 O 显著 


图 4 IBI BAS YR] T Fr BUR 328 75 ZR BICI AX ERE LEE A 


Fig.4 Linear trends of seasonal reference evapotranspiration at the meteorological stations over the Tibetan Plateau 


in different seasons during different periods 


参考 蒸 散 的 线性 变化 率 平均 值 分 别 为 10.52 mm - 
(103), 11. 29 mm + (10a) !, 7. 75 mm + 
(10a) ~' 6.84 mm + (10a) ^, FERRE 5 ZR BK 
呈 增 加 趋势 ,显著 增加 区 域 集 中 在 祁连山 (图 4b)。 
E .秋季 参考 蒸 散 大 致 在 高 原 33°N 以 南 ( 唐古拉 山 
以 南 区 域 ) 增 加 趋势 突出 ,增幅 普遍 在 10. 01 ~ 


50.00 mm + (10a) ^ Za] (K 4d 4f), BBR 
We LE FESS 1B 4 a EO A ae AS i 
小 (图 4h) ; 

总 之 ,1971 一 1996 年 高 原 参 考 蒸 散 的 减少 和 
1997—2014 年 参考 蒸 散 的 增加 主要 表现 在 春季 和 
夏季 。 
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T 旱 区 
2.2 青藏 高 原 参 考 蒸 散 变 化 趋势 的 归 因 
2.2.1 影响 参考 蒸 散 的 气候 因子 变化 “影响 青藏 


高 原 参考 燕 散 的 气候 因子 变化 趋势 存在 明显 的 区 域 
差异 (图 5)。 其 中 ,1971 一 1996 年 青藏 高 原平 均 气 
温 \ 日 照 时 数 的 线性 变化 率 为 0.15 C+ (10a) 7 

-32.79 h - (10a)”“。 这 期 间 高 原 气 温 整 体 呈 现 增 
加 趋势 (图 5a) ;日 照 时 数 下 降 的 站 点 集中 于 青藏 高 


原 东 南部 (图 5c)。 而 在 1997 一 2014 年 , 上述 要 素 
1971 一 1996 年 
80°E 90°E 100°E 
40°N 
35°N 
30N] 线性 变化 率 / 
[© + (102) J 
e <-0.50 
25^N19 029-000. eosl-100 @ 不 显著 
6001-050 6-100 Of 


40°N 
35°N 
BO°N | 线性 变化 率 


.0 
9.9 on | 50.01 ~ 100.005 
° ü 50. o 1 
25°N c 610001-150.00 © 
O001-5000 — eL150.01 


40°N 
35°N 
30°N 线性 变化 率 / 
。S-!。(10a)-!] 
d 9.30 
ooa E 
40°N 


Se 
35°N 


30?N 
PRETE (IY l 
on | 9 000. o0 
25'N16.003-0.00 @0.03 ~ 0.06 EET 
60.00-0.08 © >0.06 
80°E 90°E 100°E 


图 5 青藏 高 原 不 同时 段 气 候 因 子 的 线性 
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的 线性 变化 率 平均 值 为 0.34 "C + (10a) 7! - 4.27 

h + (10a) -'。 这 期 间 高 原 中 南部 气温 增加 显著 [ 普 
遍 高 于 0.50 *C + (10a) ^J ,高 原 北部 昆仑 山 一 阿 
尔 金 一 祁连山 附近 沿线 则 时 降温 态势 ,但 大 部 分 未 
通过 显著 性 检验 (图 5b)。 高 原 中 南部 增 温 应 该 与 
局 地 日 间 云 量 的 急速 减少 有 关 ' 。 同 期 ,高 原 33"N 
以 北 日 照 时 数 减少 趋势 显著 ,在 藏 南 谷地 附近 则 有 
所 增加 (图 5d) 。 


1997 一 2014 年 
80°E 90°E 100°E 


线性 变化 率 / 
p Ce d J 


$045 030 
0-0.29~0.00 6051-00 Ert 
2001-050 ecT” of 


° 100.00 
-99.99 ~ 50.00 50.01 ~ 100.005 
si 99 ~ 0.00 6100.01 ~ 150.00 ots 
00.01 B000" @> 150.01 


线性 变化 率 / 
(m * s?! * (102)! J 


0 
© 0.31 ~ 0.60 Tra 
e > 0.61 显著 


see Hele) 7 

$ 505 003 

$008. 000. ooos oo6 "ad 
90.00 ~ 003 @>0.06 显著 
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Fig.5 Linear trends of climatic factors over the Tibetan Plateau in different periods 
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1971—1996 年 ,青藏 高 原 各 站 点 风速 、 相 对 湿 
度 的 线性 变化 率 平均 值 分 别 为 -0.25 mss” - 
(10a) ~' 0.01 + (10a)”。 该 阶段 高 原 整 体 风 速 下 
降 显 著 ,在 高 原 北部 边缘 尤为 显著 (图 Se ) ;相对 湿 
度 显 著 增 加 区 域 则 主要 分 布 在 青藏 高 原 约 33°N 以 
南 ( 图 5g)。 突 变 后 的 1997—2014 年 ,青藏 高 原 各 
站 点 风速 .相对 湿度 的 线性 变化 率 平 均值 为 0. 07 
m*s^!-(10a) ! fl 0.04 - (10a)”。 其 中 ,青藏 
高 原 东 南部 风速 上 升 趋势 明显 (图 5f) ,这 不 同 于 突 
变 前 高 原 整 体 风速 减 小 的 态势 ;而 相对 湿度 则 在 整 
个 青藏 高 原 呈 现 全 面 下 降 趋 势 ,特别 是 在 高 原 南部 
雅鲁藏布江 中 下 游 区 域 (图 5h) ,这 与 之 前 相对 湿度 
的 增加 形成 鲜明 对 比 。 突 变 前 的 风速 减 小 应 该 是 由 


蒸 散 的 贡献 率 存 在 显著 时 空 分 异 。 突 变 前 (1971 一 
1996 年 ) ,75 个 站 点 逐步 回归 方程 R 平均 值 高 达 
0.92 ,显示 出 拟 合 优 度 较 高 。 气 候 因 子 的 标准 化 回 
归 系数 的 绝对 值 能 揭示 各 气候 因子 对 参考 蒸 散 变 化 
的 相对 重要 性 。 影 响 高 原 参 考 蒸 散 的 气温 日 照 时 
数 、 风 速 .相对 湿度 的 标准 化 回归 系数 75 个 站 点 平 
均值 分 别 为 12. 20、15.57、27.12、-15.29( 表 1)。 
可 见 , 风 速 是 影响 高 原 年 参考 蒸 散 的 关键 因子 。 风 
速 的 标准 化 回归 系数 绝对 值 居 首位 的 站 点 有 35 个 
(图 6a)。 而 风速 在 该 时 期 内 普遍 减 小 (图 Se), 
速 变 小 促使 参考 蒸 散 减少 。 据 此 ,1971 一 1996 4E, 
风速 减 小 对 高 原 参 考 蒸 散 的 急剧 下 降 起 决定 性 作 
Ho 日 照 时 数 作 用 次 之 ,其 大 范围 减少 (图 5c) 导 致 


大 气 环 流 减弱 "造成 的 ,而 纬 向 风 增 强 对 突变 后 
的 风速 变 大 有 重要 贡献 ” 。 相 对 湿度 受 水 汽 含量 
和 气温 共同 控制 ,高 原水 汽 含量 受 南亚 季风 季节 性 
进退 造成 的 水 汽 输送 变化 的 影响 "” ,而 高 原 气 温 变 
化 则 受 云 一 辐射 反馈 作用 的 控制 。 

2.2.2 气候 因子 对 参考 蒸 散 的 贡献 率 青藏 高 原 
气候 因子 (气温 \ 日 照 时 数 、 风 速 、 相 对 湿度 ) 对 参考 


高 原 参 考 蒸 散 减 少 。 相 反 ,相对 湿度 的 增加 (图 5g) 
则 使 得 高 原 参 考 蒸 散 增 加 。 

突变 后 (1997 一 2014 4E.) ,75 个 站 点 的 逐步 回 
归 方程 R^ 平均 值 高 达 0. 95 ;气温 日照 时 数 .风速 、 
相对 湿度 的 标准 化 回归 系数 75 个 站 点 平均 值 分 别 
为 13.83 14.87 .20.57 和 -20.73( 表 1)。 显 然 , 相 
对 湿度 和 风速 是 影响 该 阶段 高 原 参 考 蒸 散 变 化 格局 


表 1 影响 青藏 高 原 年 和 季节 参考 蒸 散 的 气候 因子 的 标准 化 回归 系数 


Tab.1 Standardized regression coefficients of climatic factors affecting the annual and seasonal reference 


evapotranspiration over the Tibetan Plateau 


1971—1996 年 


1997—2014 年 


年 春季 夏季 秋季 冬季 年 春季 夏季 秋季 冬季 
T 12.20 6.68 5.78 4.11 4.81 13. 83 7.53 6.28 4.84 5.14 
SD 15.57 5.47 6.30 4.47 1.34 14. 87 5.45 10.23 4.24 1.32 
U 27.12 8.42 8.61 6.61 6.29 20.57 7.13 6. 56 5.27 5.21 
RH -15.29 -6.83  -10.51 -5.44 -6.23 -20.73 -7.71 -7.42 -7.83 -5.58 


TEST 表示 气温 ;SD 表示 日 照 时 数 ;U 表示 风速 ;RH 表示 相对 湿度 。 下 同 。 


80°E 90°E 100°E 


40°N:- 
35°N 
30°N 


25°N 


(a) 1971—19964E 


(b) 1997—20144F 


6 影响 青藏 高 原 年 参考 蒸 散 的 气温 (T) .日 照 时 数 (SD) 风速 (U) 、 相 对 湿度 ( RH) 的 标准 化 回归 系数 绝对 值 


Fig.6 Absolute values of the standardized regression coefficients of annual mean temperature(T) ,annual mean relative sunshine 


duration( SD) ,annual mean 2 m height wind speed(U) , and annual mean relative humidity( RH) ,with annual 


reference evapotranspiration 
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图 7 不 同时 段 影响 青藏 高 原 各 季节 参考 蒸 散 的 气温 (T) 日照 时 数 (SD) 、 风 速 (U) 相对 湿度 (RH ) 的 标准 化 回归 系数 


Fig.7 Absolute values of the standardized regression coefficients of annual mean temperature(T) ,annual mean relative 


sunshine duration( SD) ,annual mean 2 m height wind speed(U) , and annual mean relative humidity( RH) , with annual 


reference evapotranspiration in different periods 


的 关键 气候 因子 。 考 虑 到 31 个 站 点 相对 湿度 标准 
化 回归 系数 绝对 值 最 大 (图 6b) ,再 结合 相对 湿度 在 
这 一 阶段 内 普遍 变 小 (图 5h) ,因此 ,相对 湿度 的 减 
小 对 高 原 参 考 燕 散 增加 的 贡献 率 最 大 。 此 外 ,该 阶 
段 风速 在 高 原 川 滇 藏 交界 处 增强 (图 5f) ,使 之 成 为 


高 原 参 考 蒸 散 增加 的 第 二 大 贡献 因子 。 

1997 年 前 后 2 个 阶段 驱动 高 原 参 考 蒸 散 季 节 
变化 的 气候 因子 格局 发 生 了 明显 转型 (图 7)。 
1971—1996 年 ,春秋 ,冬季 参考 蒸 散 减少 的 首要 因 
子 为 风速 的 下 降 , 春 、 秋 季风 速 的 影响 区 域 比较 相 
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似 ,大 致 在 高 原 约 33°N 以 北 区 域 ,冬季 风速 则 影响 
高 原 主体 ( 表 1, 图 7a、7c 7g)。 已 有 研究 表明 ,青藏 
高 原 高 空 南 一 北 气压 差 显 车 减弱 ,造成 我 国 大 气 环 
流 减弱 ,进而 影响 风速 降低 。 此 外 ,气温 是 影响 冬季 
参考 蒸 散 的 重要 因子 ( 表 1 ,图 7g) 。 夏 季 影 响 高 原 
参考 蒸 散 最 重要 的 影响 因子 为 相对 湿度 的 增加 ,其 
主导 区 域 为 高 原 主体 ;其 次 为 风速 的 减 小 ( 表 1, B 
7c) ,特别 是 高 原 北部 昆仑 山 一 阿尔 金山 一 祁连山 
沿线 附近 。 

1997—2014 年 ,春秋 ,冬季 参考 燕 散 增加 的 主 
导 因 子 为 相对 湿度 的 减少 , 春 、 秋 ,冬季 影响 区 域 比 
较 相 似 , 其 影响 范围 为 高 原 主体 区 域 ;同时 ,气温 上 
升 亦 促进 了 冬季 参考 燕 散 增加 ( 表 1 ,图 7b、.7f7p)。 
夏季 影响 高 原 参考 蒸 散 最 重要 的 影响 因子 为 日 照 时 
数 ( 表 1 ,图 7d) ,主要 影响 区 域 为 除了 高 原 北 部 少 
数 站 点 以 外 的 高 原 主体 区 域 。 需 要 指出 的 是 ,青藏 


高 原 上 空 夏 季 水 汽 含量 和 云 量 直 接 影响 日 照 时 
By 


3.1 结论 

基于 修正 的 FAO 56 Penman-Monteith 公式 和 青 
藏 高 原 75 个 台 站 1971—2014 年 逐日 气象 观测 资 
料 ,估算 了 高 原 40 a 来 的 参考 蒸 散 ,深入 分 析 了 参 
考 蒸 散 变化 趋势 的 转折 特征 及 发 生 转折 的 主要 原 
因 。 分 析 揭 示 了 高 原 参 考 蒸 散 在 1997 年 发 生 趋势 
转变 ,从 之 前 的 下 降 态 势 转 为 显著 上 升 趋势 。 风 速 
减 小 是 影响 1971—1996 年 高 原 参 考 蒸 散 减 小 的 主 
SEES. ,风速 的 主导 作用 遍布 整个 高 原 ,在 高 原 北部 
尤为 强烈 ;相对 湿度 下 降 则 是 导致 1997 一 2014 年 高 
原 参 考 蒸 散 增加 的 最 大 贡献 因子 。 不 同 阶 段 影响 青 
藏 高 原 参 考 蒸 散 变化 的 气候 驱动 因子 存在 显著 的 时 
空 差异 。 此 外 ,1997 年 前 后 2 个 阶段 驱动 高 原 参考 
燕 散 季节 变化 的 气候 因子 格局 发 生 了 明显 转型 。 一 
方面 ,1971 一 1996 年 ,对 春 、 秋 冬季 高 原 参 考 蒸 散 
减少 的 最 大 影响 因子 为 风速 减 小 , 且 春 、 秋 季风 速 影 
响 区 域 相似 , 均 在 高 原 约 33°N 以 北 区 域 , 冬 季风 速 
则 影响 高 原 主体 ;1997 一 2014 年 ,春秋 、 冬 季 参 考 
蒸 散 增 加 的 最 大 影响 因子 为 相对 湿度 的 下 降 ,其 影 
响 区 域 主 要 在 高 原 主体 区 域 。 另 一 方面 ,影响 高 原 
夏季 参考 蒸 散 的 主导 因子 由 1971—1996 年 相对 湿 
度 的 增加 转变 为 1997 一 2014 年 日 照 时 数 的 增加 。 


3.2 讨论 

XP HB RBS AL EAMH T 1997 年 发 
A HEE Ur ESTE DOT, 3 rr] THD ve Jie s a 2 Ze LZ 
发 量 的 时 空 变化 表现 出 较 高 的 一 致 性 ,增幅 有 所 差 
异 "”。 同 时 ,1997 年 也 是 青藏 高 原 变 暖 向 停滞 的 
转折 点 ,这 说 明了 高 原 大 气 环流 系统 的 动态 转折 
过 程 。 归 根 结 底 ,参考 疗 散 的 转变 全 面 反 映 了 青藏 
高 原 气候 系统 的 内 在 变化 。 

1997 年 之 前 ,影响 青藏 高 原 参 考 芋 散 最 关键 的 
气候 变量 是 风速 ,这 与 以 前 的 研究 结果 "一 
致 。 亚 洲 季 风 系 统 的 减弱 有 可 能 是 导致 青藏 高 原 风 
速 减弱 的 主要 原因 "” 。 但 是 ,目前 已 有 的 研究 对 其 
他 气象 变量 的 影响 作用 尚 存在 分 卜 。 笔 者 研究 认 
为 ,相对 湿度 日 照 时 数 和 气温 均 对 高 原 参 考 蒸 散 的 
减少 起 着 重要 的 作用 。 其 中 ,日 照 时 数 与 太阳 辐射 
关系 密切 ,高 原 上 太阳 辐射 降低 主要 是 由 于 水 汽 含 
量 和 深层 云层 的 增加 ,这 反 过 来 又 与 全 球 变 暖 和 对 
流 有 效 势能 的 增加 有 关 "" ,太阳 辐射 影响 高 原 东南 
部 参考 蒸 散 ,这 已 形成 广泛 共识 “  ”。 此 外 ,青藏 
高 原 日 最 低 气温 增加 显著 ,而 日 最 高 气温 增加 趋势 
不 显著 ,这 使 得 日 较 差 减 小  。 参 考 蒸 散 的 变化 
主要 受 白天 温度 影响 ,所 以 日 平均 气温 的 增加 不 能 
简单 的 说 明 是 参考 蒸 散 增加 的 原因 。 相 对 湿度 被 认 
为 影响 作用 其 微 ,然而 本 文 认为 ,相对 湿度 可 能 
是 影响 参考 燕 散 的 重要 因子 , 它 的 影响 机 理 需 要 进 
一 步 研 究 。 
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Change and Attribution of Reference Evapotranspiration over the 
Tibetan Plateau during the Period of 1971 — 2014 


WANG Bu-wei'^, ZHANG Xue-qin' 
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2. University of Chinese Academy of Sciences, Beijing 100049 , China) 


Abstract: Itis indispensable to explore the spatiotemporal variation of reference evapotranspiration over the Tibet- 
an Plateau and its driving factors so as to deepen the understanding of regional hydrological process and its response 
to climate change. The daily data from 75 meteorological stations over the plateau were utilized to analyze the spati- 
otemporal change of evapotranspiration based on the modified formula of FAO 56 Penman-Monteith, and the contri- 
bution of climatic factors to the variation during the period of 1971 —2014 was discussed. The main conclusions are 
summarized as follows: Annual reference evapotranspiration over the Tibetan Plateau tended to an increase, espe- 
cially in the region south from 33?N from 1997 to 2014 with an average linear trend of 40. 16 mm - (10a) ~'. It, 
however, decreased remarkably at the rate of - 27.07 mm - (10a) ~! from 1971 to 1996, and the dominant climat- 
ic factor resulting in the decrease was the reduction of wind speed, especially in Northern Plateau. Reduction of 
relative humidity played a crucial role in increasing the reference evapotranspiration from 1997 to 2014. In addi- 
tion, the biggest contribution factor to the seasonal variation of reference evapotranspiration shifted from the reduc- 
tion of wind speed (from 1971 to 1996) to the decrease of relative humidity ( from 1997 to 2014) and the dominant 
factor to the summer variation of reference evapotranspiration shifted from the increase of relative humidity (from 
1971 to 1996) to the increase of annual mean relative sunshine duration (from 1997 to 2014). 
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